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LATTICE CALCULATIONS OF SEMILEPTONIC FORM 

FACTORS 



J.M. Flynn 

Department of Physics and Astronomy, University of Southampton 

We review the results of lattice QCD calculations of form factors for semileptonic 
decays of D and B mesons. We also mention results for semileptonic decays of b 
baryons and the rare radiative decay B — » K**y. 

This short review deals primarily with lattice QCD results for semileptonic 
decays of D and B-mesons, but also covers the rare radiative decay B — > K*j 
as well as semileptonic decays of 6-baryons to charmed baryons- Many details 
are left out: for more information see the recent review in Rcf. El for example. 

1 Semileptonic B — ► D and B — ► D* Decays 

Semileptonic B — > D* , D decays are used to determine the CKM matrix ele- 
ment V c b- In the helicity basis the decay of a pseudoscalar meson to another 
pseudoscalar meson mediated by the vector component of the weak current is 
described by two form factors denoted / + and f°. For the decay of a pseu- 
doscalar meson to a vector meson, both vector and axial components contribute 
and there are four independent form factors, V, Aq, A\ and A2. Expressions 
for these form factors can be found in Ref.ld for example. 

Heavy quark symmetry (HQS) is rather powerful in constraining these 
heavy-to-heavy quark transitions (see Ref.B for a recent review). In the heavy 
quark limit all the form factors are described by a single universal Isgur-Wise 
function, £(u>), which contains all the non-perturbative QCD effects. Here, 
lu = v-v', the product of the initial and final meson 4-velocities. Vector current 
conservation implies the important result that the IW function is normalized 
at zero recoil, £(1) = 1. To extract \V c b\ from B — > D* decays, one extrap- 
olates the product F{uj)\V c b\ to the zero recoil point, u = 1. Here, T is the 
"physical form factor" , given by the IW function combined with perturbative 
and power corrections. One needs a theoretical evaluation of ^(1) to complete 
the extraction. In practice, one expands, 

H")=m [i-/s 2 (^-i)+c(c-i) 2 +■■■]. a) 

The slope p 2 differs from the slope p 2 of the IW function itself by heavy quark 
symmetry violating corrections □, p 2 = p 2 + (0.16 ± 0.02) + power corrections. 
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Experimental data currently show a rather wide variation in the slope and 
intercept of the extrapolation. Lattice calculations can help by providing in- 
formation on the slope of the form factor. 

To discuss lattice results for the shape of the IW function, it is convenient 
to work with a set of form factors which in the heavy quark limit either vanish 
or are equal to the IW function. These are h+(u>) and h-(w) for B — > D and 
hv(<jj), ^Ai(w), hA 2 (to) and /ia 3 (w) for B — > D*. One writes 

hi(u) = (a, + (3i{u) + 7*M)£(w) (2) 

with a+y y Ai,A 3 — 1 and a~,A 2 = 0. The /3, and ji denote, perturbative and 
power corrections (in l/mf, iC ) respectively. Luke's theorem □ states that 

7+^(1) =0(AWm» t ). ( 3 ) 

The principal difficulty for lattice calculations is to separate the physical 
heavy quark mass dependence due to power corrections from the unphysical 
one due to mass-dependent discretization errors. One must also address the 
question of lattice-to-continuum matching. We illustrate our discussion using 
the analysis procedure applied for h+ by the UKQCD collaboration B. This 
form factor is protected by Luke's theorem at zero recoil and, for degener- 
ate (Q=Q') transitions, conservation of the vector current Qj^Q provides the 
further constraints: 

P+(l;niQ,mQ) = 0, 7+(l; m Q ,m Q ) = 0. (4) 

The correct vector current normalization can be fixed by requiring the meson 
to have electric charge 1. We therefore define the continuum form factor by, 

h+{uj;mQ,mQ>) = (1 + £+(1; m Q , m Q i)) -j— r- , (5) 

where h}"{uj; mg, rriQi) is the un-normalized form factor calculated directly 
in the lattice simulation. This definition partially removes discretization er- 
rors and also removes ^-independent power corrections while maintaining the 
known normalization conditions. If the remaining power corrections are small, 
then h+(uj)/ (l + /?+(u;)) is effectively the IW function, This is convenient 

for extracting £((*>), but the definition of Eq. ^ precludes a determination of 
the zero-recoil power corrections. These corrections should be small, being 
suppressed by two powers of the heavy quark mass. However, applying an 
analogous procedure to the Ha 1 (w) form factor relevant for B — > D* decays 
will not allow the 1/mj? corrections to ^(1), one of the dominant theoretical 
uncertainties, to be determined. 
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Table 1: Slope of the IW function of a heavy meson. The table indicates which function 
from Eq. (^) has been fitted and which form factor has been used in the extraction. The 
systematic error in the UKQCD results incorporates the variation from fitting to all forms 
in Eq. ra). BSS note that , is 12% smaller than p^, but do not quote a separate result. 



Yr 



Pl,d 



LANL 
UKQCD 
UKQCD 
BSS 
BSS 



96 0.97(6) 

95 0.9(1)® 

94 0.9(f) (?) 
93 
93 



1.2(f)© 

1-2(1)® 
1.24(26)(36) 

1.41(19)(41) 



fit 

NR 

NR 

NR 

lin 

NR 



using 

h+ 
h + 

h+ 



UKQCD confirmla that their results for h + /(l + (3 + ) are indeed indepen- 
dent of the heavy quark masses and hence demonstrate, within the available 
precision, that there is an_IW function. Moreover, a similar analysis for hAi 
reveals the same function El, so the IW function appears to be universal. Ex- 
trapolating to the light (u, d) and strange (s) quark masses and fitting to 

reNR^)^^ I exp(-2(p 2 -l)^i) NR 
£(u>) = < 1 - p 2 (uj - 1) linear (6) 

[ 1 - p 2 {lo — 1) + |(cj — l) 2 quadratic 

gives lattice determinations of the slope of the IW function, as listed in Table [l| 
Since 'the' IW function is different for different light degrees of freedom, the 
results in the table are labelled with subscripts u, d or s as appropriate. 

Figure [I] shows a comparison of the UKQCD lattice resultsa with B — » D* 
data (in 1995) from CLEO A fit is made to the experimental data for 
iVyj^w). The lattice calculations cannot distinguish the u> dependence of 
T{lo) from that of £(u>) and hence p 2 u d and p 2 are not distinguished. The slope 
of the IW function is constrained to the lattice result in the fit so that the only 
free parameter is \V c b\J-(l). The result isB 

|Vc6|^(l) = 0.037+ }±!+|. (7) 

We should also mention B -=^fK semileptonic decays, which are beginning 
to be measured experimentally 113'Eil with good precision, despite the helicity 
suppression in dY{B — > DWi)/du>. The differential decay rate depends on 
both h + and h—. However, h- is rather poorly determined to date in lattice 
calculations, so that it is difficult to evaluate the 0(\/mQ) corrections. 

Direct lattice calculations of the IW function-are being undertaken using 
discretizations of the heavy quark effective theoryt 2 ], but the results are not yet 
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Figure 1: Fit of the UKQCD lattice results for |V c ;,|.7 : '(aA@ to the experimental data from 

the CLEO collaboration^]. 

useful for phenomenology. An interesting theoretical feature of this approach 
is the formulation of the HQET at non-zero velocity in Euclidean spaceES. 

2 fo-Baryon Semileptonic Decays 

There are lattice results from UKQCD for the form factors in semileptonic 
6-baryon decays to charmed baryons, in particular: — * K c e + v and —> 
S c e + ^. Heavy quark symmetry is again very predictive for these decays. 

It is convenient to use six functions of lu = v-v' as form factors for the 
weak current matrix element: ^1,2,3 for the vector part and Gi,2,3 for the axial 
vector part. At leading order in HQS they all either vanish or are given by a 
single universal baryonic Isgur-Wise function, normalised at to — 1: -Fi(a;) = 
Gi(oj) = £,(o-i), F23 = ^2,3 = 0. Luke's theorem guarantees that X^-^W an d 
Gi(l) have no 0(l/mgC)) corrections, while vector current conservation for 
Q=Q' imposes a normalization condition. This means that measuring certain 
ratios of form factors at w ^ 1 to those at ui — 1 fixes the lattice-to-continuum 
normalization and partially removes discretisation errors, as was the case for 
heavy-to-heavy meson decays described above. G\{oj) and turn out 

to be insensitive to 1/mgo corrections, so they are used to determine the 
IW function. Then G2,3 and i*i,2,3 are used to fix the l/mgo corrections, 
allowing the physical form factors to be extracted. The results for the slope of 
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Figure 2: Lattice results for D — > K and D — > K* semileptonic decay form factors at zp*p 
momentum transfer, and comparison with experimental results form the survey in Ref. c3. 
The shaded bands show our summary values, indicating which results they are based on, as 

discussed in the text. 



the baryonic IW function and for the partially integrated decay rates are: 

A b -> A c t 



P 

12 , dT 

CLLO— = 

aw 



1.21,1 



1.5: 



£=,(, — » 

;6 | 2 10 13 s 
1.6±t|K b | 2 10 13 s- 



i- 



1.4+ t|K 

5 1Kb 



7/ 



for A b 
for E b 



A c Z+t/ 
S r l+v 



(8) 
(9) 



3 Semileptonic Z? Decays 

Semileptonic Z) — » K, K* decays provide a good test for lattice calculations 
since the relevant CKM matrix element V cs is well constrained in the standard 
model. The form factors for D — > irl + vi and D — > pl + vi are also computed. 
Charm quarks are light enough to be simulated directly (though one needs to 
be wary of mass-dependent discretization errors) . Furthermore, strange quarks 
can also be simulated directly, so for D — > K, K* decays there is only one quark 
for which a chiral extrapolation needs to be performed. For semileptonic D- 
meson decays the whole physical phase space can be sampled (and beyond 
to unphysical, negative, q 2 ) 1 while keeping the spatial momenta of the ini- 
tial and final state mesons small in order to minimise momentum-dependent 
discretization errors. 

Although lattice calculations measure the q 2 dependence of the form fac- 
tors, we follow the standard practice of quoting values at q 2 — 0. In contrast 
to the case for B decays, this involves an interpolation and so is relatively well 
controlled. Lattice results for the D — > K^*' form factors are illustrated in Fig- 
ure H A full list of results for D — > and D — > n,p form factors appears 
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Table 2: Summary of lattice and experimental results for D — » K, K* and D — > it, p semitep- 
tonic decay form factors at q 2 = 0. Experimental numbers are from the survey in Ref.l 21 !. 





D K,K* 


D -f ir,p 




lattice 


expt 


lattice 


/ + (0) 


0.73(7) 


0.76(3) 


0.65(10) 


V(0) 


1.2(2) 


1.07(9) 


1.1(2) 


MO) 


0.70(7) 


0.58(3) 


0.65(7) 


A 2 (0) 


0.6(1) 


0.41(5) 


0.55(10) 



m Ref.EI These are all from quenched simulations and have not been extrapo- 
lated to the continuum. The summary in Tahle || is obtained by considering 
the more recent results from WUP0, LANL@, UKQCD0 and APE0 The 
values quoted reflect the fact that / + and A\ are the best measured while the 
D — > 7T, p form factors are smaller with slightly larger errors. 

One sees that the lattice and experimental results agree rather well. The 
lattice values for A\ and A 2 are both high compared to experiment: however, 
these depend on the correct normalization of the lattice axial vector current 
which is less well known than the vector current normalization needed for / + 
and V. In particular, the non-degeneracy of the c- and s-quarks means that 
there is no natural normalization condition to use for the weak current. This 
contrasts with the situation for heavy-to-heavy semilcptonic decays, described 
in Section [l], where one benefits from the conservation of the vector current of 
degenerate quarks. 

4 Semileptonic B — > p and B — > tt Decays and B — > K*j 

The heavy-to-light semileptonic decays B — > p and B rp>|7r are now being used 
experimentally to determine the V u b matrix element LilE3. Several-groups, hapj 
evaluated form factors for these decays using lattice simulations li3'li§E3G3 c3 
(see the recent review in Ref.c3). We will also consider the rare radiative decay 
B — > K*j which is related by heavy quark and light flavour symmetries to the 
B — > p semileptonic decay. The physical B — > K*"f amplitude is determined 
by the value of form factors Ti(0) or T 2 (Q) at the on-shell point q 2 =0. With 
the definitions used here (see, for example, Ref. El), Ti(0) = ^2(0). 

Heavy quark symmetry is less predictive for heavy-to-light decays than for 
heavy-to-heavy ones. In particular, there is no normalization condition at zero 
recoil corresponding to the condition £(1) = 1, so useful in the extraction of 
V c b- This puts a premium on results from nonperturbative techniques, such as 
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lattice QCD. HQS does, however, give useful scaling laws for the form factors 
as the mass of the heavy quark varies at fixed uj. Moreover, the heavy quark 
spin symmetry relates the B — > V matrix elements E3c3 (where V is a light 
vector particle) of the weak current and magnetic moment operators, thereby 
relating £>° — > p + l~vi and B — > K*j, up to SU(3) flavour symmetry breaking 
effects. For fixed u> the scaling laws for the form factors given by HQS are as 
follows: 

Z(9 2 M) I fixcd e = M Vf v + + + (10) 

where / labels the form factor, M is the mass of the heavy-light meson and 
is a calculable leadiag logarithmic correction. The leading M dependences, 
M v f , are listed in Ref.u for example. These relations can be used to extrapolate 
lattice calculations with quark masses around the charm mass to the B mass. 
In the limit M — > oo we also have the relations 

A 1 (q 2 (oj)) = 2iT 2 (q 2 (cu)), V{q 2 {u)) = 2T 1 (< 7 2 (u,)) I (11) 
at fixed oj. The UKQCD collaboration have checked the validity of the relations 

mH 



in Eq. ( |ll| ) finding that they are well satisfied in the infinite mass limit. 
Finally, there are also kinematic constraints on the form factors at q 2 = 0: 

/+(0) = /°(0), Ti(0)=iT 2 (0), A (0)=A 3 (0). (12) 

To control discretization errors in lattice simulations we require that the 
three-momenta of the B, 7r and p mesons be small in lattice units and therefore 
we determine the form factors only at large values of momentum transfer q 2 . 
Experiments can already reconstruct exclusive semileptonic b — > u decays (see, 
for example, the review in Ref. E3) and in the future we can expect to compare 
the lattice form factor calculations directly with experimental data at large q 2 . 
A proposal in this direction was made by UKQCDcJ for B° —* p + l~i?i decays. 
They parametrize the differential decay rate distribution near q^ax by: 

where a and b are parameters, and X(q 2 ) = {m 2 B + m 2 — q 2 ) 2 — Am^m 2 . The 
constant a plays the role of the IW function evaluated at u) = 1 for heavy-to- 
heavy transitions, but in this case there ia_no symmetry to determine its value 
at leading order in HQS. UKQCD obtainN 

a = 4.6±8;i±0.6GeV, b = (-8+ 1) x 1(T 2 GeV 2 . (14) 
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q 2 (GeV 2 ) 

Figure 3: Differencial decay rate as a function of q 2 for the semileptonic decay B° — > o^l~vi. 
taken from Ref. Ej. Point s a re measured lattice data, solid curve is fit from Eq. (hsh with 
parameters given in Eq. (tfj). The dashed curves show the variation from the statistical 
errors in the fit parameters. The vertical dotted line marks the charm endpoint. 

The result for a incorporates a systematic error dominated by the uncertainty 
ascribed to discretization errors and would lead to an extraction of \V u b\ with 
less than 10% statistical error and about 12% systematic error from the theoret- 
ical input. The prediction for the dT/dq 2 distribution based on these numbers 
is presented in Figure [|. 

We would also like to know the full q 2 dependence of the form factors, 
which involves a large extrapolation from the high q 2 values where lattice 
calculations produce results. In particular the radiative decay B — > K *7 occurs 
at q 2 = 0, so that existing lattice simulations cannot make a direct calculation 
of the necessary form factors. 

An interesting approach to this extrapolation problem has been applied 
by Lellouchca for B° — » ir + l~Di. Using dispersion relations constrained by 
UKQCD lattice results at large values of q 2 and kinematical constraints at 
q 2 = 0, one can tighten the bomids on form factors at all values of q 2 . The 
results (at 50% CL — see Ref. for details) are 

/ + (0) = 0.10-0.57, T{B° -nr + rui) = 4.4-13 |V^|ps -1 . (15) 

In principle, this method can also be applied to B — > p decays. Recently, 
BecirevicEl has applied the method for B — > K*j, using APE lattice results 
as constraints. These dispersive analyses can provide model-independent re- 
sults, but, unfortunately, the resulting bounds are not very restrictive when 
constrained by existing lattice data. 

For now we must rely on model input to guide q 2 extrapolations. Wc 
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Table 3: Lattice results for B — » n,p semileptonic decays. Rates in units of |V u (,| 2 ps 







Yr 


B° — > 
Rate 


/ + (0) 


B°-^p+ 
Rate 


I Vl 

V(0) 


M°) 


A 2 (0) 


UKQCD 
WUP 
APE 
ELC 


2E 
!( 

L7 


97 
97 
95 
94 


8.5(?i; 

8±4 
9±6 


0.27(11) 
0.43(19) 
0.35(8) 
0.30(14) (5) 


16.5(H) 

12 ±6 
14 ±12 


0.35(g) 
0.65(15) 
0.53(31) 
0.37(11) 


0.27© 
0.28(3) 
0.24(12) 
0.22(5) 


0.26(|) 
0.46(23) 
0.27(80) 
0.49(21)(5) 



can ensure that any model is consistent with HQS, as shown in Eq. (|l0|), to- 
gether with the kinematic relations of Eq. (|l2|). Even with these constraints, 
however, current lattice data do not by themselves distinguish a preferred <rL: 
dependence. More guidance is available from light-cone sum rule analyses til 
which lead to scaling laws for the form factors at fixed (low) q 2 rather than 
at fixed u> as in Eq. (|Tol). In particular all form factors have a leading M~ 3 / 2 
dependence at q 2 =0. It is important to use ansatze for the form factors com- 
patible with as many of the known constraints as possible. 

Lattice results for B° =i tt + 1~Di, -BfL— > p + l~vi and B — > K*-f are reported 
in Tables | and |. ELC Ell and APE E3 fit lattice data for the semileptonic 
decays at a single value of q 2 to a simple pole form with the appropriate pole 
mass also determined by their data. For the f° and A\ form factors, this 
is consistent with heavy quark symmetry requirements, kinematic relations 
and light-cone scaling relations at q 2 — 0_but for the other form factors it is 
not simultaneously consistent. The WUPcUl results are found by scaling form 
factors at q 2 = from results with quark masses around the charm mass to 
the 6-quark mass. However, the scaling laws used do not follow the light-cone 
scaling relations. There are also preliminary results for heavy-to-light form 
factors from FNAL, JLQCD and a Hiroshima-KEK group (see the reviews in 
Refs.E3u3) and the different lattice calculations are in agreement for the form 
factors at large q 2 where they axe measured. 

The latest UKQCD study l3 uses models consistent with all constraints, 
including the light-cone sum rule scaling relations. In lattice calculations one 
has the freedom to adjust hadron masses by tuning the quark masses used in 
the simulation. UKQCD use this freedom to perform a combined fit for all 
the B — > V form factors (where V denotes a light final state vector meson) 
simultaneously, first with V = p and then with V = K* . They obtain form 
factors for B° — > p + l~Di from the first fit and for B — > K*j from the second. 
The combined fit in the K* case is illustrated in Figure [|, which demonstrates 
the large extrapolation needed to reach q 2 — 0. 
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i igure 4: UKQCD0 fit to the lattice predictions for Aq, Ai, V, Ti and T2 for a K* meson 
final state assuming a pole form for A\. The dashed vertical line indicates g^ax' 

Table 4: Lattice results for B -> IC7. Values for T(0) = Ti(0) = iT 2 (0) are quoted only 
from models satisfying the light-cone sum rule scaling relation at q 2 = 0. 







Yr 


T(0) 




UKQCD 


25 


97 


0.16(f) 


0.25(2) 


LANL 


33 


96 


0.09(1) 




APE 


34 


96 


0.09(1)(1) 




BHS 


35 


94 


0.101(10)(28) 


0.325(33)(65) 



Our preferred resultsJbr B° —* ir + l~vi and B° — > p + l~Di come from the 
UKQCD constrained fitsEj. Their values for the form factors extrapolated to 
q 2 = agree well with light-cone sum rule calculations, which work best at low 
q 2 . The fitted form factors also agree with experimental results for the rates 
and ratio-of-rates of these semileptonic decays. However, we emphasise that 
the extrapolated form factors are no longer model independent. 

Table | lists the values of T(0) = Ti(0) = £T 2 (0) for B -> K*j, together 
with the directly measured T2{q 2 aax ). All groups find that T2 has much less 
q 2 dependence than T\. The table lists results from form factor fits satisfying 
the light conc_sum rule scaling relation at q 2 = 0. Our preference is to quote 
the UKQCD @ result, T(Q) = 0.16(f) (with statistical error only), from the 
combined fit described above. Using this one cari evaluate the ratio (at leading 
order in QCD and up to 0(l/m 2 ) corrections @) R K * = T(B -> K*^)/Y(b -> 
S7) = _Lfi(3)%. This is consistent with the experimental result 18(7)% from 

cleoeI. 
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